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Abstract 


We have discussed the physical mechanisms through which changes In the 
boundary forcings of SST, soil moisture, albedo, sea Ice, and snow Influence 
the atmospheric circulation. The slowly changing boundary forcings can 
increase the predictability of monthly means because their effects on quasi - 
stationary flow patterns and statistics of synoptic scale disturbances appears 
to be potentially predictable. Changes In the boundary forcings produce changes 
In the moisture sources and dlabatic heat sources which In turn change the 
atmospheric circulation. The magnitude and the structure of the atmospheric 
response due to changes in any boundary forcing depends upon the existence of a 
suitable large scale flow which can transform the boundary forcing into a three 
dimensional heat source, which in turn can change the large scale flow and its 
stability properties. The structure of the large scale flow also affects the 
propagation characteristics of the influence which determines whether the effect 
is local or away from the source. 

We have presented results of numerical experiments conducted with the 6LAS 
climate model to determine the sensitivity of the model atmosphere to changes 
in boundary conditions of SST, soil moisture, and albedo over limited regions. 

It is found that changes in SST and soil moisture in the tropics produce large 
changes in the atmospheric circulation and rainfall over the tropics as well as 
over mid-latitudes. Although the area occupied by the land surfaces is small 
compared to the ocean surfaces, the fluctuations of soil moisture can be very 
important because the dlabatic heat sources have their maxima over the land, 
and therefore, even small fluctuations of soil moisture can produce large changes 
in the total diabatic heating field. The natural variability due to day to day 
weather fluctuations is very large in the middle latitudes, and therefore, changes 
in the mid-latitude atmospheric circulation due to changes in the boundary 


forcings at middle and high latitudes have to be quite large to be significant. 
It Is suggested that large scale persistent anomalies of SST, snow and sea Ice, 
under favorable conditions of large scale flow, can produce significant changes 
In the mid-latitude atmospheric circulation. It Is also likely that time 
averaged mid-latitude circulation can have additional predictability due to the 
Influence of tropical boundary forcings. 

We have also presented observational evidence to show that Interannual 
variability of atmospheric fluctuations Is significantly different from the 
Intra-annual variability, and therefore, we conclude that part of the inter- 
annual variability Is due to the Influence of boundary forcings. Since the 
tropical spectra Is found to be redder compared to the mid-latitude spectra, 
the tropical flows may be potentially more predictable. 
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1. INTRODUCTION 

In a previous paper (Shukla, 1981) we examined the predictability of the 
initial conditions without external forcings, and It was shown that the dynamical 
predictability of the observed planetary wave configurations Is sufficiently 
long that the predicted monthly means are significantly different from the 
monthly means due to random perturbations in the initial conditions. It was 
further suggested that In addition to the dynamical predictability of the 
monthly means there can be additional predictability due to the influence of 
the boundary forcings. This paper examines the Influence of slowly varying 
boundary forcings at the earth's surface in determining the monthly mean 
circulation anomalies In the atmosphere. Figures la and lb show the schematic 
distribution of sea surface temperature (SST), soil moisture, surface albedo, 
snow cover and sea Ice. There Is sufficient observational evidence to assume 
that the rate of change of anomalies of these boundary forcings is slower than 
the corresponding atmospheric anomalies, and therefore for a limited period 
of time (viz. a month or sometimes even a season) these can be considered as 
"external" forcing to the atmospheric circulation. In the present paper, we 
will confine our discussion only to these five boundary forcings. A discussion 
of the influence of truly external forcing due to solar variability is beyond 
the scope of the present paper. Similarly, effects of slow changes in the 
composition (viz. CO 2 ) and optical properties of the atmosphere wi il not be 
considered because they appear to be too slow to affect monthly prediction. 

One straightforward approach to determining the influence of boundary condi- 
tions on monthly mean prediction could be to carry out actual forecast experi- 
ments using the observed initial and boundary conditions. To our knowledge, no 
such forecast experiment has yet been carried out in which the global distribu - 
tion of observed boundary forcings was used. It is only recently that global 
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observing systems capable of determining global boundary conditions have been 
established, and still It Is a formidable (but feasible) task to determine 
a global distribution of accurate boundary forcings. Another approach to 
determine the Influence of the boundary forcings Is to carry out a series of 
controlled sensitivity experiments with realistic general circulation models 
and to determine the regional distributions of time averaged response. This 
has been done using a variety of models and a variety of anomalous boundary 
forcings. A review of all such experiments carried out by many investigators 
Is not possible 1r this paper. Here we summarize the results of only a few 
selected numerical experiments carried out with the Goddard Laboratory for 
Atmospheric Sciences (GLAS) climate model. 

It should be pointed out that most of the numerical experiments carried 
out In the past were not especially designed to Investigate the question of 
monthly predictability and, therefore, there was no detailed examination of 
the response during the first 30 days. The main Interest was the difference in 
the simulated mean climate after the model has equilibrated to the altered 
boundary forcings. Such experiments are also quite useful for understanding the 
basic mechanisms and establish the role of boundary forcings in the interannual 
variability of monthly and seasonal means. In earlier studies, the choice of 
the geographical location of the anomaly (of SST or sea ice, etc.) was mostly 
determined by an analysis of the past observations, and naturally the past 
observations were analysed only for those areas for which data were conveniently 
available. On the basis of such observational studies and subsequent numerical 
experiments, several preliminary conclusions have been drawn about the possible 
relationships between the boundary forcing anomalies and the atmospheric cir- 
culation. These studies do not rule out the possibility that SST anomalies 
over other areas (hereto unexamined) will not be equally Important. In fact. 
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It Is quite likely that all the boundary forcing, If considered globally, will 
produce more realistic and systematic response in the atmospheric circulation. 

One of the main contentions of this paper Is to emphasize that further sensitivity 
and predictability studies using observations, simple linear models and complex 
GCMs should also be carried out for global distributions of boundary forcing 
anomalies. 

Before presenting the results of actual numerical experiments, we will 
briefly mention some mechanistic considerations which suggest a physical basis 
for the influence of the boundary forcings. 

1) Changes in the boundary forcings directly Influence the location and 
Intensity of the diabatic heat sources which drive the atmospheric circulation. 

The forcing at the boundary itself is generally not sufficient to produce 
significant changes in the atmospheric circulation; however, under favorable 
conditions of large scale convergence and divergence, the boundary effects get 
transmitted to the interior of the atmosphere and a thermal boundary forcing 
gets transformed into a three-dimensional heat source which can be quite effec- 
tive in influencing the dynamical circulation. The effectiveness of a boundary 
forcing in changing the atmospheric circulation therefore strongly depends 

upon its ability to produce a deep heat source and the ability of this influence 
to propagate away from the source. Since both of these factors are determined 
by the structure of the large scale dynamical circulation itself, the response 
of a given boundary forcing can be very different depending upon its size and 
geographical location, and upon the structure of the large scale circulation. 

2) Boundary forcings of SST, soil moisture, surface albedo, snow and ice 
not only affect the heat sources and sinks, but they also affect the sources 
and sinks of moisture, which in turn affect the latent heat sources. 

3) Existence of nonlinear multiple equilibrium states for a prescribed 
external forcing suggests that even weaker anomalies of boundary forcings. 


under favorable conditions, can produce significant anomalies In atmospheric 
circulation, and therefore actual response may be stronger compared to the one 
estimated from linear theories. 
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4) It Is known that the Inconsistency between the observed Initial condi- 
tions and the prescribed stationary forcings due to mountains and dlabatlc 
heat sources can manifest Itself as erroneous propagating transient components 
(Lambert and Merllees, 1978; Shukla and Llndzen, 1981). It Is therefore desirable 
that fo" actual dynamical prediction from observed initial conditions, the observed 
global distribution of boundary forcings be used correctly. This will reduce 
the Inconsistency between the Initial conditions and the forcing, and therefore 
can reduce the growth of error of prediction. 

In Part I we examined predictability for prescribed nontluctuatlng boundary 
forcings, and the predictability was determined by error growth rate and error 
saturation value. We can Identify the following reasons due to which changes 
In boundary forcings can Influence these classical predictability parameters. 

a) Changes In the boundary forcings can change the intensity and geograph- 
ical location of synoptic scale instabilities. For example, changes in the 
amplitudes and phases of planetary waves can affect the storm tracks, which can 
affect the error growth and predictability In a particular region. 

b) Changes In the boundary forcings can alter the saturation value of the 
error. As described earlier, the saturation value of the error depends upon 
the equilibration mechanisms which are different for different circulation 
regimes and also for different scales of Instabilities. Boundary forcings can 
alter the equilibration level of the dominant fluctuations, resulting in a 
significant change in the time averaged mean circulation. It should be pointed 
out that the effects can be very different for different parameters. For 
example. If for a given value of boundary forcing, the amplitude of the wave 
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disturbances Is much less than Its amplitude for another value of the boundary 
forcing, and If heavy rain or snow falls only during half of the life of each 
wave disturbance, and If the rate of rainfall Is proportional to the intensity 
of the disturbance, then although the time averaged value of pressure or temper- 
ature will not be very different for two values of boundary forcing, the rainfall 
will be very different. This suggests that sometimes It may be quite useful 
to predict only the variance of a particular parameter. 

2. SENSITIVITY OF MODEL ATMOSPHERE TO CHANGES IN BOUNDARY FORCINGS 

In this section, we have summarized the results of several sensitivity 
studies carried out with global general circulation models to determine the 
Influence of prescribed changes in SST, soil moisture, surface albedo, snow and 
sea Ice, etc. 

2.1 Sea Surface Temperature 

The factors which determine the Influence of sea surface temperature (SST) 
anomalies can be briefly summarized as follows. 

(I ) The magnitude and the spatial and temporal structure of the anomaly . 
Considering the linear response of the atmospheric system to diabatlc forcing, 
the magnitude and the spatial scale of the anomaly can directly affect the 
response to the atmosphere. Large, persistent anomalies can produce a larger 
response than small fluctuating anomalies. The magnitude of the ar..,naly is 
also Important In determining the nonlinear Increase of evaporation and sensible 
heat fluxes. 

(II) Normal sea surface temperature. Due to the nonlinearity of the Clausius- 
Clapeyron equation, a 1° positive anomaly over a normal SST of 30° produces a 
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much larger change In saturation vapor pressure than the same 1° anomaly super- 
imposed upon a normal temperature of 20°. It Is partly for this reason that 
anomalies In low latitudes can produce larger responses compared to similar 
anomalies in middle latitudes. The final response also depends upon the back- 
ground SST field. 

(Ill) The latitude of the anomaly. Because of the smallness of the Coriolis 
parameter, horizontal temperature gradients In the tropics produce a much larger 
thermal wind than In middle latitudes. Due to the lack of geostrophic constraint 
In low latitudes, thermal anomalies produce much larger convergence than in 
middle latitudes (Hoskins and Karoly, 1981; Webster, 1981). 

(1v) Circulation regime . The potential response of a given anomaly strongly 
depends upon the structure and dynamics of the circulation regime In which the 
anomaly Is embedded. For example, a warm anomaly In the areas of large scale 
convergence (viz. ascending branches of Hadley and Walker cells) will be more 
effective than a comparable anomaly In the area of divergence. Similarly, In 
the middle latitudes the effect of a SST anomaly will strongly depend upon the 
location of the anomaly with respect to the phase of the prevailing planetary 
wave configurations. 

(v) Instability mechanism . The time required for the atmosphere to feel 
the effect of the SST anomaly also depend: upon the most dominant Instability 
mechanism which determines the generation of a deep heat source due to surface 
anomaly. In tropical latitudes where CISK Is the primary driving mechanism, a 
conditionally unstable atmosphere may respond rather quickly to a warm SST 
anomaly, whereas In mid-latitudes where the primary driving mechanism Is barocllnic 
Instability a given SST anomaly would affect the vertical shear and, therefore, 
the growth rates of barocllnlcally unstable waves. 

(vl) Structure of zonal flow . Once a SST anomaly has produced a heat 
source, the structure of the prevailing zonal flow is of crucial importance n 
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determining the propagation characteristics of the disturbances produced by the 
heat source. Tropical Influences can affect the mid-latitude circulation by 
Rossby wave propagation or by changing the Intensity of the Hadley cell and mid- 
latitude zonal flows which can Interact with the mid-latitude thermal and orographic 
forcings. 

In summary, the Influence of boundary anomalies on the atmospheric circulation 
depends upon the existence of a favorable dynamical environment In which the 
surface forcing can be transformed Into a three-dimensional heat source and the 
ability of this Influence to propogate away from heat source. A warm anomaly 
In the tropics enhances evaporation and Increases the moisture flux convergence 
which Is the main contributor to the enhanced precipitation over the anomaly. 
Increased evaporation lowers the lifting condensation level, Increases the buoyancy 
of the moist air, accelerates the deep convective activity, and Increases the 
latent heating of the atmospheric column which further reduces the sjr'ace 
pressure and enhances moisture convergence. In the tropics, SST anomalies can 
also produce considerable effects away from the anomaly by modifying the areas 
of convergence and divergence. For example, if the intertroplcal convergence 
zone remains stationary over a very warm SST anomaly for a considerable length 
of time, those areas where ITCZ would have moved in Its normal seasonal march 
will experience severe droughts. Similarly, a warm SST anomaly can change the 
location and Intensity of the Walker circulation and enhanced ascending motion 
associated with the warm SST anomaly can produce reduced precipitation in the 
adjoining areas. This suggests that the effects of SST anomalies can be very 
nonlinear for particular regions, depending upon the location of the region 
with respect tc the ascending branches of Hadley and Walker circulations. This 
nonlinearity can disappear for averages over very large area*. SST anomalies 
can also affect precipitation over distant areas by altering the moisture 
supply for the region. 
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The Influence of SST anomalies for mid-latitudes Is different than that 
for the tropics. SST anomalies In mid-latitudes, If they can produce a deep heat 
source, can change the quasi -stationary wave patterns which In turn can affect 
the location and Intensity of the storm tracks. Therefore, the mid-latitude 
SST anomalies have a cons' Arable potential to produce distant effects. This 
potential, however, 1$ not fully realized for several reasons. The normal 
ocean temperature In mid-latitudes Is relatively cold; the moist convection Is 
not well organized and efficient enough to produce aeep heat sources; due to 
strong geostrophlc balance even large gradients In temperature do not produce 
large convergence, and finally, the natural variability of the mid-latitude 
atmosphere Is so large that small effects due to boundary forcings cannot he 
distinguished from synoptic weather fluctuations. 

Most of the observational aid numerical studies that have been carried out 
so far have only looked Into the influence of regional SST anomalies. We do 
not have enough observational and numerical experimental results to describe the 
effects of hemispheric and global scale SST anomalies. We hope that future 
studies will examine the effects of very large scale SST anomalies. Kere we 
shall present the results of a few sensitivity studies carried out by global 
general circulation models to determine the Influence of regional SST anomalies. 

2.1.1. Effect of Arabian Sea sea surface temperature anomaly on Indian monsoon 
raWall 

Several observational studies have suggested that SST anomalies over the 
Arabian Sea can be one of the Important boundary forcings which determine the 
monsoon rainfall and therefore affects Its Interannual variability. A numerical 
experiment was carried out by Shukla (1975) using the Geophycical Fluid Dynamics 
Laboratory (6FDL) model to test the validity of observed correlations. A cold 
SST anomaly shown In Figure 2a (anomaly run) was Imposed over the cllmato’ofp' ’ 
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SST (control run) and the model was Integrated for both cases. Figure 2b shows 
the model simulated rainfall of the Indian region for anomaly and control Inte- 
gral on«. It is found that due to cold SST anomalies over the Arabian Sea 
monsoon rainfall is reduced over India. A similar experiment was carried out 
by Washington et al . , (1977) using the National Center for Atmospheric Research 
(NCAR) model; however, their results were contrary to both the numerical results 
of the 6FDL model and observed correlations. We have examined the results of 
the Goddard Laboratory for Atmospheric Sciences (GLAS) climate model in which 
a warm SST anomaly very similar in pattern to the anomaly shown in Figure 2a 
was used. Figures 3a and 3b show the model simulated rainfall for anomaly and 
control runs averaged over two areas shown in Figure 3c. The GLAS model was 
integrated for three different initial conditions, but very similar SST fields 
over the Arabian Sea. The curves labeled ANOMALY in Figures 3a and 3b show 
the average rainfall for three anomaly runs and the vertical bars denote the 
standard deviation among the three runs. It can be seen that warm SST anomalies 
over the Arabian Sea, as prescribed in these experiments, can produce enhanced 
monsoon rainfall over India. We do not know, however, whether such large 
anomalies are actually observed. 

The apparent disagreement between the results of Washington et al . , (1977) 
using the NCAR model, and our results using the GFDL and GLAS models can be 
explained by examining the low level monsoon flow as simulated by the three 
models and shown in Figure 3d. The low level monsoon flow as simulated by GFDL 
and GLAS models are more realistic compared to the NCAR model. In the NCAR 
model simulation, the air parcels flowing over the Arabian Sea hardly reach 
the Indian region as there is a strong but unrealistic southward flow before 
the monsoon current reaches the Indian coast. This could be one of the possible 
reasons why the NCAR r :.al did not show significant response over India. This 
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example Illustrates the Importance of a realistic simulation of the mean climate 
by the dynamical model used for the sensitivity studies. In order to be able 
to detect the effects of changes In the boundary forcings, and for the possible 
use of such models for prediction of time averages, the model should be able to 
simulate the mean climate accurately. 

2.1.2. Effect of tropical Atlantic sea surface anomalies on drought over 
northeast Brazil 

Moura and Shukla (1981) examined the monthly mean SST anomalies over the 
tropical Atlantic durlnj March and rainfall anomalies over northeast Brazil 
during March, April and May. They found that the most severe drought events 
were associated with the simultaneous occurrences of warm SST anomalies over the 
north tropical Atlantic and cold SST anomalies over the south tropical Atlantic. 
They also carried out numerical experiments to test the sensitivity of the 6LAS 
climate model to prescribed SST anomalies over the tropical Atlantic. It was found 
that the SST anomaly patterns, which resemble the observed ones during drought 
years, produced an intensified convergence zone (ITCZ), enhanced rainfall and 
low level cyclonic circulatK.i to the north, and reduced rainfall and anticy- 
clonlc circulation to the south. 

Figure 4a shows the SST anomaly used for the numerical experiment; the 
magnitude was chosen to be comparable to the maximum values observed during 
the 25 year period (1948-72). Figure 4b shows the 15 day running mean time 
series of dally rainfall averaged over the areas A and B (shown in Figure 4a) 
for control and anomaly runs. Area A Includes the region of warm SST anomaly 
and area B contains the northeast Brazil region and the neighboring oceans with 
the cold SST anomaly. The rainfall over area A increases due to the warm SST 
anomaly and a shift of the ITCZ occurs from area B in the control run to area 
A In the anomaly run. Although the difference in rainfall is not systematic 
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after 60 days, for days 20-60 the anomaly run has consistently less rainfall 
than the control run over northeast Brazil region. Figure 4c shows the differ- 
ence of the first 60 day mean meridional circulation averaged between 50°W and 
5°E. The anomalous meridional circulation shows an ascending branch with 
maximum vertical motions between 5° and 10°N and a descending branch to the 
south of the equator. A thermally direct local circulation is established with 
its ascending branch at about 10°N and its descending branch over northeast 
Brazil and adjoining oceanic regions. The driving for the anomalous circula- 
tion is provided by convection and latent heating associated with warmer SST 
anomalies over the northern tropical Atlantic, and cooling associated with 
colder SST anomalies in the southern tropical Atlantic. The combined effects 
of thermally forced subsidence and the reduced evaporation and moisture flux 
convergence produces severe drought conditions over northeast Brazil. It should 
be pointed out that descending motion influences a larger region to the east 
and west of northeast Brazil, however, the effect is most seriously felt over 
northeast Brazil. It was also noticed that the anomalous meridional circula- 
tions exhibited significant changes in the middle latitudes of both hemisphere. 

Figure 4d shows the autocorrelation function for SST anomalies at 15°S, 

5°W and at 15°N, 45°W. The strong persistence of SST anomalies in these two 
areas suggests the potential for prediction of drought over northeast Brazil. 

2.1.3 Effect of Equatorial Pacific SST Anomaly on Tropical and Extratropical 
Ci rcui at ions . 

Horel and Wallace (1981) and Rasmusson and Carpenter (1982) have presented 
observational evidences of remarkable relationships between SST anomalies in 
tropical Pacific and a variety of atmospheric fluctuations including the Southern 
Oscillation and the Northern Hemispheric middle latitude circulations. "Warm 
episodes of equatorial Pacific SST anomalies are associated with the negative 


phase of the Southern Oscillation, weakening of easterlies In equatorial central 
Pacific, enhanced precipitation at equatorial stations east of 160°E, an Intensi- 
fied Hadley cell In Pacific sector, and a deepening and southward displacement 
of the Aleutian low" (Horel and Wallace, 1981). Analyses by Rasmusson and 
Carpenter (1982) have shown that during the month of December, a small warm SST 
anomaly appears along the Peru coast which rapidly Increases to its peak value 
along the Peru coast during March and April of the following year. In the month 
of February of the following year, the warm SST anomaly disappears rapidly along 
the Peru coast, but the warmest SST anomalies are observed in equatorial central 
Pacific. It is these warm SST anomalies over the central Pacific, which occur 
over a climatologically warmer ocean surface which Is also an area of large 
scale convergence, that produce a deep tropical heat source whose effects can 
propagate to the extratroplcal latitudes. Therefore, although the signal for a 
warm SST anomaly during February of a given year could be traced back to warm 
SST anomalies over the Peru coast fourteen months earlier, it is the enhanced 
heating associated with the SST anomaly during February that produces significant 
changes in the winter circulation of the northern mid-latitudes. Hoskins and 
Karoly (1981) used a multi-level linear primitive equation model to show that 
dlabatlc heat sources in the tropical regions can produce significant stationary 
responses in the middle latitudes If the zonal flow is favorable for the propa- 
gation of Rossby waves. 

Shukla and Wallace (1983) have conducted sensitivity experiments with the 
GLAS climate model to study the response of SST anomalies in equatorial Pacific. 
Figure 5a shows the average of SST anomalies observed during the months of 
November, December and January of 1957-58, 1965-66, 1969-70, and 1972-73, 
provided to us by Dr. Rasmusson rf NOAA. The GLAS climate model was integrated 
with (anomaly run) and without (control run) the SST anomalies shown in Figure 5a. 
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Figure 5b shows the difference (anomaly minus control) field for 300 mb geopo- 
tential height averaged for days 11-25. It Is seen that a difference of about 
300 meters over North America and about -90 to -150 meters over Pacific and 
Atlantic is in agreement with results of observational studies as well as the 
results of linear models. Upper level highs to the north of the anomaly and a 
series of lows and highs further north and west are manifestations of Rossby 
waves propagating from the heat source. An examination of day-to-day changes 
from day 1 through 30 showed that this particular pattern was already established 
during days 5-10. However, it takes several days before model physics can 
generate a deep heat source above the warm SST anomaly. It should be pointed 
out that although the results of the linear models and this general circulation 
model are similar, there is a basic difference for their applicability to 
dynamical prediction of monthly means. The simple models prescribe the diabatic 
heating fields, whereas the general circulation model generates a diabatic 
heating field due to the presence of the warm SST anomaly, and therefore for a 
model to be useful as a prediction tool, its physical parameterizations mist 
be able to transform the surface boundary forcing into a deep heat source. As 
mentioned earlier, the propagation of the tropical influences depends upon the 
structure of the zonal flow (influences cannot propagate across the zero wind 
line). In the present experiment, although the SST anomaly was centered right 
over the equator, the zonal flow was favorable for the tropical effects to 
propagate to middle latitides. 

We have also examined the changes in the model simulated Hadley and Walker 
cells due to the prescribed SST anomalies. Figure 5c shows the difference 
between the Hadley cell for the anomaly and the control run averaged for days 
6-25. It is seen that the zonal 1y averaged Hadley cell intensifies due to 
warm SST anomalies over the equatorial Pacific. The Hadley cell over the Pacific 
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sector alone Is Intensified even more. It was also noticed that In association 
with stronger Hadley cells, the westerly zonal flow was also stronger between 
20° and 30°N. Figure 5d shows the difference (anomaly minus control) for the 
model s mulated Walker cells averaged between 6°N and 6°S for days 6-25. 

Anomalous ascending motion occurs between the longitude sector 160°E - 160°W 
and descending motion occurs between the longitudes 140°E - 160°E. This is 
consistent with the observational evidence of enhanced precipitation east of 
160°E being associated with warm SST anomalies. 

A remarkable aspect of the results of this experiment Is that the pres- 
cribed SST anomaly was an algebraic mean of four different episodes and there- 
fore the prescribed maximum warm anomaly at any grid point was perhaps less 
than observed during any Individual year. However, it produced significant 
changes within the first 30 days. The results also suggest that the SST ano- 
malies can be of Importance even for medium range (5-15 days) forecasting. 

2.1.4 Effect of north Pacific SST anomaly on the circulation over North America . 

As It was pointed out earlier, due to the large day-to-day variability in 
the middle latitudes. It is difficult to detect the influence of SST anomalies. 
However, if a large scale SST anomaly of large magnitude persists for a long 
time, it can also produce significant effects in the midlatitude circulation. 
During the fall and winter of 1976-77, SST In the north Pacific was characterized 
by abnormally cold temperatures In the central and western portions of the 
northern Pacific with a warm pool located off the west coast of the U.S. (Figure 
6a). Namlas (1978) has suggested that the northern Pacific SST anomalies may 
have been one of the multiple causes of the abnormally cold temperatures in the 
eastern North America during the 1976-77 winter. We have carried out numerical 
experiments with the 6LAS climate model to test this hypothesis. It should be 
noted that although the pattern of the SST anomalies observed during January 
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1977 was similar to the one shown In Figure 6a, the magnitude of the anomaly 
was only slightly more than half of the magnitude shown In the figure. It 
should also be pointed out that during January 1977 a warm SST anomaly over the 
equatorial Pacific was also observed and it is likely that the circulation of 
North America could have been affected by combined effects of both anomalies. 

In this experiment, we have attempted to determine the effect of north Pacific 
SST anomaly only. We have integrated the GLAS climate model with six different 
initial conditions and climatological mean SST (control run) and for two of 
these Initial conditions we have also integrated the model with imposed SST 
anomalies (anomaly run). The results that we show here are the differences 
between an average of six control runs and two anomaly runs. 

Figure 6b shows the difference (anomaly minus control) for the 700 mb 
temperature averaged for days 15-45. As expected, negative temperature ano- 
malies are found over the maximum negative SST anomaly and a warm 700 mb temp- 
erature anomaly is found over the warm SST anomaly. The most remarkable 
feature of this map, however, is the occurrence of two other centers of nega- 
tive 700 mb temperature anomaly, one along 75°W and the other along 70°E. The 
cold 700 mb temperature anomaly over northeast North America is of largest 
magnitude. This feature suggests a planetary wave response due to north Pacific 
SST anomaly. The ratio between the temperature differences shown in Figure 6b 
and the standard deviation among the six control runs with climatological SST 
shows (Figure 6c) that the colder temperatures at 700 mb over northeast North 
America are as significant as the ones over the cold SST anomaly. The difference 
map (anomaly minus control) for the 500 mb geopotential height field shows (Figure 
6d) that colder temperatures over northeast North America were due to anomalous 
northerly flow of cold air which in turn was caused by anomalous positive 


16 


geopotential height anomalies to the west and negative geopotentlal height 
anomalies to the east. There are discrepancies between the simulated and the 
observed anomaly during the winter of 1977, because the response Is one quarter 
wavelength out of phase. However, a detailed examination of the day-to-day 
evolution of the flow for the control and anomaly runs revealed that a persistent 
blocking type of flow configuration existed In the anomaly run for more than 
15 days (Chen and Shukla, 1983). Such persistent blocks did not occur In the 
control case. We are not quite sure If the blocking event generated In the 
anomaly run was due to presence of the SST anomaly. However, on the basis of 
the results of these experiments and the results of several observational 
studies. It is our conclusion that for favorable structures of the large scale 
flow, large scale SST anomalies In the mid-latitudes can be important In deter- 
mining the anomalies of the mid-latitude circulation. 

2.2 Soil Moisture 

The annual average relnfall for the global continents is estimated to be 
about 764 mm of which 35-40% (266 mm) runs off to the oceans (Baumgartner and 
Reichel , 1975). Assuming no secular trends In the annual mean global soil 
moisture, this suggests that the annual and global mean evaporation from tlie 
land surfaces alone Is more than 60% of the annual and global mean precipitation 
over the land. The percentage Is even higher during the local summer for 
several regions. This suggests that the evaporation from the land surfaces is 
a very Important component of the global water budget and hydrological cycle. 
However, It does not necessarily follow that the water evaporated from the 
land Is Important In determining the rainfall over the land. For example, all 
the water evporated from the land could be advected away to the oceans before 
It recondenses and rains. In that case. It will affect the moisture budget 
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and evaporation only over the oceans, which In turn will, of course, affect 
the moisture supply for rainfall over the land. In order that the evaporation 
from the land affects the rainfall over the land, it Is necessary that the 
prevailing dynamical circulation be such that the land evaporated moisture 
recondenses and falls as rain before being advected away. That will depend upon 
the geographical location of the region under consideration, the prevailing 
advectlve velocity, the structure and Intensity of the convergence field, and 
the vertical distribution of moist static energy which determines the nature 
of moist convection. 

The role of soil moisture is twofold. First, it determines the rate of 
evaporation and, therefore, the moisture supply, and second, it influences the 
heating of the ground which determines the sensible heat flux and affects the 
dynamical circulation by generation or dissipation of heat lows. The interac- 
tion between the heat lows, generated by solar heating of the ground In the 
absence of soil moisture, and associated circulation and rainfall Is further 
complicated by the fact that the maintenance and the intensification of the 
low pressure areas is largely Influenced by the latent heat of condensation. 

For example, if the soil is saturated with water, and the evaporation is equal 
to the potential evapotranspi ration, there will be maximum possible supply of 
moisture to the atmosphere. Whether Increased evaporation will also increase 
the rainfall will depend upon the structure of dynamical circulation and 
prevailing flow patterns. If the rate at which the moist static energy (cpT + 
gz + Lq) is advected away from the region is larger than its accumulation rate. 

It will not lead to any Increase In the rainfall. For the other extreme situa- 
tion, when the soil is completely dry, and there is no evaporation from the land, 
there may be a reduction In the rainfall due to reduced evaporation. However, 
if the heating of the land produces Intense low pressure areas which can converge 
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moisture from the surrounding oceans, the rainfall may not necessarily decrease, 
and If the convergence of moisture Is large enough rainfall may even Increase. 
The mechanism will cease to operate, however, once the rain starts falling 
because the soil will not be dry anymore. 

Since the net dlabatlc heating of a vertical atmospheric column Is maximum 
over the tropical land masses (Figures 7a and 7b), It Is quite likely that 
small fractional changes In these tropical asymmetric heat sources could produce 
considerable changes In the planetary scale circulations of the tropical as 
well as the extra-tropical atmosphere. Therefore, In spite of relatively 
smaller earth surface area being covered by land, soil moisture effects could 
be as Important as SST anomaly effects. It should be noted, however, that the 
soil moisture effects strongly depend upon the season and latitude because 
during the winter season In high latitudes, solar radiation reaching the ground 
Is not large enough to be Important for surface energy budget. We summarize 
here the results of two numerical experiments, carried out by Shukla and Mlntz 
(1982) which have demonstrated that evaporation from land can significantly 
affect the rainfall over land. 

2.2.1. Influence of global dry-soil and wet-soil on atmospheric circulation 

We have carried out two 60 day integrations of the GLAS climate model: in 

one case, there Is no evaporation from the land surface ("dry-soil" case) and 
In the other case the evaporation from the land Is equal to the model calculated 
potential evapotransplratlon (“wet-soil" case). These two cases are qualita- 
tively similar to no vegetation and completely vegetated earth surface. For 
simplicity of Interpretation of the results, albedo of the soil was not altered 
for the two experiments. 

Figure 7c shows the global maps of mean July rainfall difference (dry soil- 
wet soil). Over most of the continental regions, with the exception of tht 
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major monsoon regions, July rainfall for dry-soil case has decreased (by about 
40-50*) compared to the wet-soil case. This shows that the evaporation from 
the land Is an Important component of the rainfall over the land. The exception 
over India Is even more Interesting because the solar heating of land for dry- 
soil produced such Intense low pressure and convergence that the loss of moisture 
from land evaporation was more than compensated by the Increased moisture flux 
convergence from the neighboring oceanic region. The Increased moisture flux con- 
vergence leads to Increased heating of the vertical air column due to the latent 
heat of condensation which maintains and Intensifies the surface low. This 
happened prominently for the unique monsoonal circulation over India for which 
oceanic moisture was brought In from the ocean by the monsoon current. This 
could also occur In several other regions with monsoonal flow patterns but 
perhaps the coarse grid resolution of the model could not resolve the local 
Influence over other regions. 

Figure 7d shows the ground temperature difference (dry soil-wet soil) for 
July. Ground temperatures for dry-soil are warmer by more than 20°-30°C. For 
the dry-soil case, most of the radiation energy goes to heat the ground and to 
Increase the sensible heat flux, whereas, for the wet-soil case, most of the radia 
tlon energy goes to evaporate the water which can later release the latent heat of 
condensation and heat the vertical air-column. 

Figure 7e shows the surface pressure difference (dry soil -wet soil) for 

# 

July. There are Intense low pressure areas over the continents for dry-soil 
and the mass removed from over the land Is found to produce high pressure areas 
over the oceans. The location and the Intensity of the high pressure areas 
over the oceans are determined by the nature of dynamical circulations, the 
most Important of which are the mid-latitude stationary wave response to a 
highly anomalous dlabatlc forcing, and the modified Hadley and Walker circulations 


20 


due to changes In the Intensity and location of the tropical heat sources. 

We have also calculated the natural variability (not shown) of the GLAS 
model and It should suffice to say that the changes shown here for dry-soil and 
wet-soil are far too large to be confused with the model variability due to 
Internal dynamics. 

Since the spatially and temporally averaged rainfall over any region 
depends upon evaporation (which depends upon soil moisture), and vertical 
distribution of the moisture flux convergence (which depends upon the nature 
of dynamical circulation, which. In turn Is determined by the direct heating 
of the ground and vertical distribution of latent heat of condensation), It 
1$ not possible to establish simple universal relationships between soil moisture, 
atmospheric circulation and rainfall. Such relationships strongly depend upon 
the dynamical circulation In the region under consideration. If the moisture 
over land Is rapidly advected away to oceans. It Is unlikely that the local 
evaporation from land will be an Important contributor to the rainfall. On 
the other hand, the Intensity and life cycle of a tropical disturbance which 
moves over land can strongly depend upon the wetness of the ground. It Is 
therefore necessary to utilize realistic physical models of the earth-atmosphere 
system to determine the sensitivity of climate to fluctuations of soil moisture. 

The experiments reported In this study are extreme examples to highlight the 
maximum bounds of Impact and to show that the Importance of soil moisture and 
vegetation for rainfall over any particular area should be calculated by 
Including the combined effects of dynamical circulation and orop-aphlc and 

i 

oceanic forcing. It has also been reported that soil moisture boundary forcing 
can be an Important factor for medium range (5-15 day) forecasting of rainfall 
and circulation. (See Mlntz (1982) for a review of several numerical experiments.) 
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2.3 Surface Albedo. 

Charney et al. (1977) had suggested that changes In the surface albedo can 
produce significant changes In the local rainfall and atmospheric circulation. 

They pointed out that this effect can be especially Important In the desert 
margin regions of subtropics. Charney et al . used an earlier version of the 
GLAS model In which the surface albedo was increased by 30% for the Sahel, the 
Thar Desert (India), and the western Great Plains of the U.S. They found that 
precipitation over the albedo anomaly regions was reduced by 10-25% within the 
first 30 days. An Increase In the surface albedo reduced the solar radiation 
reaching the ground, which In turn reduced evaporation from the ground. These 
reductions In evaporation and cloudiness increased the solar radiation reaching 
the ground, thus partly compensating for reduction In solar radiation cue to 
Increased albedo. However, since reduction In the c’oudlness also caused 
reduction In the long wave radiation emitted back to the surface from the 
cloud base, there was a net reduction In the total radiative energy coming to 
the ground. This caused a net reduction in evaporation, cloudiness and precipi- 
tation. Since these subtropical regions were not affected by large advectlve 
effects, the local changes In radiative and latent heating were accompanied by 
dynamical circulations which produced descending motion over the albedo anomaly 
regions. As It was pointed out by Charney et al. (1977), the net effect on the 
atmospheric circulation due to changes in surface albedo depends upon the relative 
magnitude of the time scale for advecting the meist static energy away from the 
region and the time scale for Its generation by evaporation and convergence. 

If the net effect of the change in surface albedo is to reduce the sensible and 
latent heating of the air. It will either decrease the low level convergence 
and ascending motion or Increase the low level divergence and descending motion. 
The change In the vertical velocity field In either case will dry the middle 
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troposphere and further reduce precipitation. These effects may not operate 
In areas of large sr.ale moisture convergence or strong dynamical Instabilities; 
however* these experiments suggested that It Is quite Important to give a 
realistic prescription of surface albedo for prediction of monthly and seasonal 
atmospheric anomalies. 

These experiments have been recently repeated by Sud and Fennessy (1982) 
using the present version of the GLAS climate model which has better parameter- 
Izatlons for evaporation and sensible heat fluxes. Sud and Fennessy have found 
that the results of experiments with the present model support all the conclu- 
sions of the earlier studies by Charney et al. An Increase In albedo produced 
a systematic decrease In rainfall over Sahel and Thar Desert. The only exception 
occurred for western Great Plains where the change was not large enough to be 
distinguished from the day-to-day variability. 

2.3.1. Effect of change In albedo over northeast Brazl 1 . 

Ir. addition to the areas chosen by Charney et al. , Sud and Fennessy also 
examined the effect of change In surface albedo over northeast Brazil. Between 
4°S and 24°$, and 32°M and 47.5°W, surface albedo for 10 model grid polnis was 
Increased from about 9% to 30%. Figure 7f shows 5 day averages of model simu- 
lated total rainfall over northeast Brazil for the anomaly (Increased albedo) and 
the control run. It Is seen that the rainfall for the first 20 days decreases 
due to an Increase In albedo. For a F day period from day 20-25 rainfall In 
the anomaly case Is more than the control. But again, from day 25-45 rainfall 
for the anomaly case Is considerably less than that for the control case. Since 
rainfall over northeast Brazil Is sometimes affected by the penetration of Southern 
Hemisphere mld-litltude disturbances, the changes during day 20-25 could be 
attributed to extratroplcal variability. 
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2.4 Snow Cover 

About a century ago, Blanford (1884) observed that, "The excessive winter 
and spring snowfall in the Himalayas is prejudicial to the subsequent monsoon 
rainfall in India." These observations were later substantiated by Walker 
(1910). A persistent anomaly of snow cover can affect the meridional tempera- 
ture gradient and, therefore, vertical shear of the large scale flow. The 
monsoon circulation is characterized by reversal of the normal temperature 
gradient between the equator and 30°N, i.e. during the Asiatic summer monsoon 
season, the equator is colder than northern India. An excessive snowfall during 
the previous winter and spring season can delay the build-up of the monsoonal 
temperature gradients because most of the solar energy will be utilized for 
evaporating the snow or for evaporating the soil moisture due to excessive 
snow. A weaker meridional temperature gradient between equator and 30°N can 
give rise to a delayed and weaker monsoon circulation. No systematic numerical 
experiments have been carried out with global GCMs to determine the physical 
mechanisms that can affect the atmospheric circulation associated with excessive 
snowfall. One reason for lack of such numerical experiments is perhaps the 
requirement of a rather long time integration (from one season to the other) 
of a global general circulation model. It requires a realistic treatment of 
albedo as well as surface hydrology because a large and deep snow cover during 
winter and spring can keep the soil wet for a longer time in the coming summer, 
and this effect must be treated accurately. It is only recently that general 
circulation models have shown some success in simulating the seasonal cycle of 
the atmosphere, and it is hoped that more systematic studies of various snow 
cover feedback mechanisms will be carried out in the coming years. 

Hahn and Shukla (1976) found an apparent relationship between Eurasian snow 
cover and Indian monsoon rainfall (shown in Figure 8). Large and persistent 


winter snow cover anomalies over Eurasia can produce a colder mid-latitude 
troposphere In the following spring which can strengthen the upper level anti- 
cyclone, slow Its northward movement over India, and give rise to delayed and 
weaker monsoon rainfall. Some recent studies by Yeh et al. , (1981) have suggested 
that Eurasian snow cover anomalies can Influence the interannual atmospheric 
variability of China. 

Namias (1962, 1978) has proposed that positive feedback can occur between 
excessive srow cover on the east coast of North America and quasi -stationary 
circulations which can be favorable for producing more snowfall. Since the 
natural variability of the mid-latitude atmosphere is quite large, it is difficult 
to determine the contribution of snow cover anomalies in producing mid-latitude 
atmospheric anomalies. The possible physical processes which have the potential 
to influence the atmospheric circulation due to snow cover anomalies can be 
briefly mentioned: 

(i) An increase in snow cover increases the albedo, and therefore reduces 
the incoming solar radiation. If there are no other feedbacks, this will 
produce colder temperatures, and therefore snow cover anomalies will tend to 
persist for a longer time. This is consistent with the observations of Wiesnet 
and Matson (1976) who showed that December snow cover for the Northern Hemisphere 
is a very good predictor of snow cover for the following January through March. 

(ii) Excessive snow cover anomalies in the mid-latitudes can produce anomalous 
diabatic heat sources which in turn can produce anomalous stationary wave 
patterns which can alter storm tracks and their frequency. 

(iii) Persistent snow cover anomalies can change the components of the heat 
balance of the earth's surface. Even after the snow has melted completely, wet 
soil can maintain colder surface temperatures for longer periods of time. 

(iv) Persistent snow anomalies can produce anomalous meridional temperature 
gradients and anomalous vertical wind shears. Snow is also a good Insulator, 
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and in the presence of deep snow cover, nighttime Inversions tend to be much 
stronger. 

2.5 Sea Ice 

Although sea ice fluctuations cover a very small fraction of the earth's 
surface, they can produce important changes in the atmospheric circulation over 
the polar regions, and for favorable dynamical structure of the large scale flow, 
these effects could also propagate to middle and subtropical latitudes. The 
local effect of sea ice anomalies is very large, because it directly changes the 
heat and moisture supply to the atmosphere. There have been several observa- 
tional studies to determine possible relationships between sea ice extent and 
atmospheric circulation anomalies (for a review of these studies, see Walsh 
(1978), and Walsh and Johnson (1978, 1979)). In a numerical experiment conducted 
by Herman and Johnson (1979) using the GLAS climate model, it was found that the 
sea ice anomalies in the Arctic regions not only affected the local circulation 
but they also produced significant differences in middle and subtropical lati- 
tudes. More experimentation with simple models and realistic GCMs are needed 
to understand the sea ice effects. Since sea ice anomalies, like SST and snow, 
also change slowly compared to the atmospheric anomalies they can be prescribed 
from observations. 

3. PREDICTABILITY OF THE TROPICAL ATMOSPHERE 

The tropical atmosphere is potentially more predictable than the mid-latitudes 
because its planetary scale circulations are dominated by monsoon circulations 
which are Intrinsically more stable than mid-latitude Rossby regime. Interac- 
tion of the large scale overturnings with the tropical disturbances (easterly 
waves, depressions, cyclones, etc.) is not strong enough to detract from the 
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predictability of the planetary scale circulations. Tropical disturbances are 
Initiated by barotroplc-barocllnlc Instabilities but their main energy source 
Is the latent heat of condensation. Although their growth rate Is fast and they 
are deterministically less predictable, their amplitude equilibration is also 
quite rapid and with the exception of hurricanes, they attain only moderate 
Intensity. The Intensity and geographical locations of the planetary scale 
circulations are primarily determined by the boundary conditions and not by 
synoptic scale disturbances. It Is reasonable to assume that frequency and 
tracks of depressions and easterly waves are primarily determined by the location 
and intensity of the planetary scales and distribution of SST and soil moisture 
fields. It is highly unlikely that the synoptic scale tropical disturbances, 
through their Interaction with planetary scale disturbances, will drastically 
alter the character of large scale tropical circulation. This is in marked 
contrast to the case of mid-latitudes where interaction between synoptic scale 
instabilities and planetary scale circulations is sufficiently strong so that 
baroclinically unstable disturbances can make the large scales unpredictable: 
the mid-latitude circulation consists of a continuum of scales whereas tropical 
circulation has a clear scale separation in terms of the frequency and the 
zonal wavenumber. 

Charney and Shukla (1981) have suggested that since large scale monsoon 
circulation is stable with respect to dynamic instabilities, and since boundary 
conditions exert significant influence on the time averaged monsoon flow, the 
monsoon circulation is potentially more predictable than the middle latitude 
circulation. This suggestion was made by examining the variability among the 
monthly mean (July) circulation of four model runs for which the boundary 
conditions were kept Identical, but the initial conditions were randomly per- 
turbed. It was found that although the observed and the model variabilities 
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were comparible for middle and high latitudes, the variability among the four 
model runs for the monsoon region was far less than the observed Interannual 
variability of the atmosphere as a whole. This led to the suggestion that part 
of the remaining variability could be due to the boundary forcings. 

We have extended the work of Charney and Shukla (1981) and compared the 
model variability for climatological and observed SST anomalies. We have 
carried out 45 day integrations of the GLAS climate model for seven different 
sets of Initial conditions and boundary conditions. In four of chese integrations, 
the climatological global SST was used. In the remaining three, the observed 
SST for 1973, 1974 and 1975 was used between 0-30°N. Figure 9 shows the plots 
of zonally averaged values of standard deviations op, og and o 0 , and ratios 
o 0 /op and o 0 /ob as a function of latitude. o p is the model standard 
deviation among predictability integrations (climatological SST and random 
perturbation In initial conditions), og Is the model standard deviation for 
SST anomalies, and o 0 is the standard deviation for 10 years of observations. 

It is seen that, in agreement with the results of Charney and Shukla, the ratio 
oo/°p is more than two in the topics and close to one in the middle latitudes. 

The new result of this study is that og (variability due to changes in SST 
boundary conditions) lies nearly halfway between o 0 and ap. This suggests 
that nearly half of the potentially predictable variability is accounted for by 
changes in SST between 0-30°N. 

These conclusions are further supported by a more comprehensive study by 
Manabe and Hahn (1981), who integrated the GFDL spectral climate model for 15 
years with prescribed but seasonally varying boundary conditions of SST, Figure 
10a (reproduced from Manabe and Hahn (1981)) shows the values of zonally aver- 
aged standard deviation of 1000 mb geopotential height for a 15 year model run 
(o m ) and observations (o 0 ). We have calculated the ratio (o 0 /o m ) from the 
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two curves of Manabe and Hahn and the ratio Is also shown In Figure 10a. It Is 
again seen that the ratio (o 0 /o ffl ) Is about two In the near equatorial regions 
and reduces to about one In the middle and high latitudes. Figure 10b (which 
Is reproduced from Figure 5.10 of Manabe and Hahn) shows the latitude and 
height cross sections of Interannual variability of geopotentlal height. It 
Is seen that the ratio between the observed variability and the model variability 
is more than three In the tropical upper troposphere. It should be pointed 
out, however, that although SST was not varying from one year to the other for 
the 15 year model run, the soil moisture and snow cover were still variable 
during different years, and therefore. It Is likely that part of the simulated 
model variability In the tropics, and perhaps even In the middle latitudes, 
could be due to the Interannual variability of soil moisture and snow cover. 

It Is reasonable to conclude that although for short and medium range the 
tropical atmosphere Is less predictable, the time averages (monthly and seasonal 
means) for the tropics have more potential predictability. Since there Is 
sufficient evidence that tropical heat sources can also Influence the middle 
latitude circulation. It Is likely that, even for mid-latitudes, the monthly 
means could be potentially predictable due to their interaction with low lati- 
tudes. 

4. POTENTIAL PREDICTABILITY OF MONTHLY MEANS AS DEDUCED FROM OBSERVATIONS 

As pointed out in the Introduction, the Interannual variability of monthly 
means Is determined by complex Interactions between the Internal dynamics and 
surface boundary forcings. It Is therefore difficult to determine their separate 
contributions by analysis of observed data. One possible approach Is to carry 
out controlled numerical experiments described earlier, provided that all the 
boundary forcings can be reasonably prescribed. 
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Hadden (1976) has examined the predictability of monthly mean sea level 
pressure over the Northern Hemisphere by comparing the variances of the observed 
monthly means with the natural variability of the monthly means. His study 
concluded that there are several areas over the Northern Hemisphere for which 
the ratio of the above two variances Is more than one and therefore there Is 
potential predictability. However, In Madden's study, the predictable signal 
was assumed to be only that part of the variance which Is above white noise 
for more than 96 day time scales. He further assumed that the boundary 
forcings do not contribute to variability for time scales shorter than 96 days. 
Since It Is well known that changes In boundary forcings can Influence atmospheric 
fluctuations within 96 days, Madden's results should be considered only as the 
lower bound of the estimate of potential predictability of the monthly means of 
the atmospheric circulations (Shukla, 1983). Moreover, Madden's study does not 
address the question of dynamical prediction from an initial state. 

Shukla and Gutzler (1983) have carried out an analysis of variance to 
compare the interannual variability of 500 mb geopotential height among 15 
(1963-77) January months, and variability within each January. It can be hypo- 
thesized that if variability among the Januaries of different years Is signifi- 
cantly larger than that due to the day to day variability within the individual 
Januaries, the excess variability could be due to boundary forcings. Since 
daily values within the Individual Januaries are not independent of each other, 
we first calculated the time Interval between Independent samples which were 
used to calculate the effective sample size and effective degrees of freedom, 
and then calculated the natural variability of January means. 

Figure 11a shows the ratio of the observed variability among the Januaries 
of different years and the natural variability. It is seen that the ratio 
between the two variances Is larger than two (In some places, even larger than 
three) over a substantial part of the Northern Hemisphere. This suggests that 
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the boundary forcings play an Important role In determining the Interannual 
variability of monthly means. As mentioned earlier, these observational studies 
also suggest that for time averages (monthly and seasonal means) the tropical 
atmosphere Is potentially more predictable than the mid-latitude atmosphere. 

This suggestion Is further supported by a comparison of autocorrelation functions 
for the tropics and mid-latitudes. Figure lib shows the zonal ly averaged values 
of autocorrelation calculated from 53 years (1925-77) of monthly mean sea level 
pressure data over the Northern Hemisphere analyzed by Trenberth and Paollno 
(1981). It Is seen that the autocorrelation drops off much more sharply for 
mid-latitudes than In the tropics. The tropical spectra are redder than the 
mid-latitude spectra and therefore potentially more predictable. 

5. CONCLUSIONS 

We have briefly reviewed the physical mechanisms through which boundary 
forcings can influence the variability and predictability of monthly means. We 
have used a realistic global general circulation model to determine the sensi- 
tivity of the model circulation to changes In the boundary conditions of SST, 
soil moisture and surface albedo. From these studies It can be concluded that 
the slowly varying boundary forcings have the potential to increase the predic- 
tability of time averages. It Is further suggested that correct specification 
of global boundary conditions of SST, soil moisture, surface albedo, sea ice, 
and snow Is necessary for successful dynamical prediction of monthly means. 

Anomalies of slowly varying boundary forcings produce anomalous sources of 
heat and moisture which In turn produce significant anomalies of atmospheric 
circulation. These effects can be either local In the vicinity of the boundary 
anomaly or away from the source If the Intervening environment is favorable for 
the propagation of local Influences. For example, tropical heat sources can 
produce significant changes In the extra-tropical circulation. 
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It should be pointed out that the anomalies of SST or soil moisture are 
neither necessary nor sufficient to produce changes in the mid-latitude circul- 
ation. For example, a warm SST anomaly In tropics can produce a local heat source, 
but If It occurs within the band of easterlies, Its Influence cannot propagate 
beyond the zero wind line. Conversely, even In the absence of SST and soil 
moisture anomalies. Internal dynamics can produce persistent anomalies of preci- 
pitation and dlabatlc heating In the tropics which, under favorable zonal flow, 
can affect the mid-latitudes. 

The main purpose of the two papers (Shukla (1981) and this paper) was to 
examine and establish a physical basis for dynamical prediction of monthly 
means. The essential reqlurements for establishing such a basis are to show 
that: a) fluctuations of monthly means are larger than can be expected due to 

sampling of day-to-day weather changes; b) there are low frequency planetary 
scale components of the circulation which remain predictable beyond the limits 
of synoptic scale predictability; and, c) Influences of the slowly varying 
boundary conditions of SST, soil moisture, snow, sea ice, etc., are large enough 
to produce significant and detectable changes in the monthly mean circulations. 

In these two papers we have presented observational and numerical evidence 
which support the above requirements. 

It is hoped that an accurate description of global boundary forcings can 
be obtained on a near-real time basis so that feasibility of experimental pre- 
diction of monthly means can be investigated using dynamical models. 
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